Dumbbell transport and deflection in a spatially periodic potential 
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We present theoretical results on the deterministic and stochastic motion of a dumbbell carried by 
a uniform flow through a three-dimensional spatially periodic potential. Depending on parameters 
like the flow velocity, there are two different kinds of movement: transport along a potential valley 
and stair-like motion oblique to the potential trenches. The crossover between these two regimes, 
as well as the deflection angle, depends on the size of the dumbbell. Moreover, thermal fluctuations 
cause a resonance-like variation in the deflection angle as a function of the dumbbell extension. 

PACS numbers: 47.61.-k, 05.40.-a, 05.60.-k 



Introduction.- The exploration of particle motion in 
microfluidic devices and energy landscapes has attracted 
considerable attention recently, in particular due to its 
great potential for applications in chemistry and biotech- 
nology [l|, [H, @, 0, @, B, 0] • With the successful genera- 
tion of optically induced potentials by holographic laser 
tweezers [H, [H , a powerful tool has been introduced for 
investigations of the classical transport of colloidal parti- 
cles through those two- or three-dimensional landscapes. 
Since the interaction between the colloids and the po- 
tential depends on the size and shape of the particles, 
such potentials may also be used in microfluidic devices 
for a spatial decomposition of different particle species. 
For instance, spherical particles, carried by a flow with 
low Reynolds number, can be deflected into a direction 
oblique to the flow lines, while transversing a spatially pe- 
riodic energy landscape [U, Q . Since the deflection angle 
is a function of the size and refraction index of the par- 
ticles this behavior may be used for particle sorting. For 
extended and dcformable objects like vesicles, instead of 
spherical particles, a similar cross-streamline deflection 
has been found in Poiseuille flow even in the absence of 
a potential 

In this context, interesting questions arise: How are 
dcformable objects, like dumbbells that are carried by a 
flow, deflected by a spatially periodic potential? What ef- 
fects are related to the interplay between the wavelength 
of the potential and the size of the dumbbell [T3|? What 
is the role of the hydrodynamic interaction between the 
two connected spherical particles? 

Our numerical investigation of the flow driven dumb- 
bell motion in spatially periodic potentials can serve as 
a simple model for other non-spherical objects, such as 
dimers , pom-pom polymers [l2j , or two small beads 
connected by polymers. We find in this system two flow 
velocity thresholds separating three different regimes of 
motion: a locked stage, transport along a potential trench 
and a stair-like motion. In addition we detect a remark- 
able dependence of the dumbbell deflection angle on the 
ratio A/6, where A is the wavelength of the periodic po- 
tential and b is the length of the dumbbell. Furthermore, 
in the presence of stochastic thermal forces, interesting 
deflection resonances have been found as a function of 
A/6. 



Model- We investigate the dynamics of a dumbbell 
that consists of two beads connected by a linear spring, 
and driven by the uniform flow u = u(cosa, sin a, 0) 
through a three-dimensional spatially periodic potential. 
The over-damped motion of the dumbbell is described by 
the Langevin equation, 



Ti = u + H 



(i,i = 1,2) (1) 



for the two bead positions ri^. The spring force, Ff = 
— V$(ri — Tz), between the beads is derived from the 
harmonic potential, 
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with the spring constant k and the equilibrium bead dis- 
tance b. The force Ff = — W(r.;) is caused by the 
periodic potential, 



V(Ti 



-V [cos (qxi) + cos (qyi) + cos (qzi)] , (3) 



with wavelength A = 2ir/q in all three spatial directions, 
and amplitude Vo- The third contribution on the right 
hand side of Eq. ([3]) is included in order to keep the 
dumbbell axis in the xy plane. In experiments such po- 
tentials can be realized by laser tweezers [1, Q, where 
Vo may be changed by varying the intensity of the laser 
beam. 

Eq. |T]) is a nonlinear function of the bead distance 
rtj due to the mobility matrix H.y, which describes the 
hydrodynamic interactions (HI) between the two beads. 
Within the Rotne-Prager approximation 13[ it has, for 
i 7^ j, the form 
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with Vij = Vi — Yj and = |ry|. The diagonal compo- 
nent of the mobility tensor, = 4 I, with unity matrix 
I, is inversely proportional to the Stokes friction coeffi- 
cient, C = 67r?7a, of a sphere with effective hydrodynamic 
radius a within a solvent of viscosity r\. 

The uncorrelated stochastic force Ff in Eq. ([1]) has 
a vanishing mean value. The amplitude of Ff is deter- 
mined by the fluctuation dissipation theorem [l4j and 
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includes the thermal energy ksT: 
<Ff(f)> 
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Eq. JTJ) has been solved numerically by an Euler method, 
as approximate analytical solutions can only be obtained 
in some special cases. In our simulations we have fixed 
the values of A = 1 and Vq = 1, while all other lengths 
and energies are given in terms of them. Furthermore, 
we choose a radius of a = 0.2 for the beads, a spring 
constant of k = 10 and a Stokes friction coefficient of 
£ = 10. Angles are given in radians with respect to the x 
axis. The time scale is set by the relaxation time of the 
spring t = C/k, therefore the velocities are given in units 
of A/t. Since the flow direction u is restricted to the 
xy plane, the periodic variation of the potential in the z 
direction keeps the dumbbell axis in a single minimum 
with respect to z, i.e. z = 0. At the initial state of the 
calculations, we align the dumbbell parallel to the flow 
lines. 

A typical situation investigated in this work is sketched 
in Fig. [T] The dotted lines in the figure indicate the 
valleys of the periodic potential that are parallel to the x 
axis (trenches parallel to the y axis are not shown). The 
streamlines of the flow field u cross these valleys with an 
angle a < it /A. This flow imposes a drag force on the 
beads of the dumbbell, which is proportional to u, and 
which depends on the bead distance r\2- If F^ is large 
enough the dumbbell moves with the mean velocity v. 
Its deflection angle is denoted by (3. 
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FIG. 1: A sample deterministic trajectory of the center of 
mass of a dumbbell (solid line) in the xy plane, which is driven 
by the flow u through a three-dimensional periodic potential 
V. The mean velocity of the dumbbell is denoted by v and 
the horizontal dotted lines indicate the valleys of V parallel 
to the x axis. 

Results.- We first consider the dumbbell motion in the 
limit of vanishing thermal noise. For small flow strengths 
u the dumbbell is captured by the potential. If the flow 
velocity u = u(l,0,0) exceeds a certain threshold value 
u c i , the dumbbell surmounts the saddles of the potential 



at a height Vq and moves along a wavy potential trench 
with a periodically changing velocity parallel to the x 
axis. The threshold velocity u c i is a function of the bead 
distance b as shown in Fig. [5] for one parameter set. 




FIG. 2: The threshold velocity u c i for driving the dumbbell 
along the x direction is shown as a function of the equilibrium 
bead distance b (solid line). The dotted line indicates the 
analytical result calculated from Eq. {7}. 

The dumbbell aligned with its axis parallel to the x 
direction and with a fixed bead distance b experiences a 
potential force related to Eq. (J3J), 



F v {x, x + b) = -— (V(x) + V(x + b)) 
ox 



(7) 



This force vanishes at bead distances b = (2n + l)A/2, 
and accordingly the required flow velocity for driving 
the dumbbell continously through the potential also van- 
ishes. However, at finite values of the spring constant k, 
corresponding to flexible objects like polymers etc., the 
dumbbell is stretched and compressed during its motion 
through the spatially modulated potential. With this de- 
gree of freedom the threshold flow velocity, u c \, still has 
pronounced minima at the distances b = (2n+ l)A/2 but 
it does not vanish as in the case of a fixed bead distance, 
especially when the HI is taken into account. This can be 
seen in Fig. [21 where the dotted line indicates the thresh- 
old velocity of a rigid dumbbell. Moreover the HI reduces 
the drag force on the two beads. Since this reduction is 
proportional to 1/t , 12 , the maxima of the curve u c i(b) 
decrease with increasing values of b. 

If the flow velocity, u = u(cos a, sin a, 0), is increased 
further, one reaches, for finite angles a and beyond a 
second threshold velocity u C 2 > u c i, another regime of 
dumbbell motion. If the y component of the drag force 
F 1 " is large enough, the beads frequently jump to a neigh- 
boring trench of the potential parallel to the x axis. Con- 
sequently the dumbbell performs a stair-like motion in 
the xy plane as indicated by the solid line in Fig. [TJ The 
mean direction v of this motion has a finite deflection 
angle (3, which depends on u as shown in Fig. [3J In the 
limit of large values of u, (3 approaches the inclination 
angle a of the flow u as expected. 
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FIG. 3: The deflection angle j3 of the dumbbell is shown as a 
function of the flow velocity u at a — 0.2 for different values 
of the bead distance 6. In the regime of f3 — the dumbbell 
moves along a potential trench. 



Since the drag force imposed by the flow on the dumb- 
bell depends on the distance between the two beads. ri2, 
the deflection angle is sensitive to the equilibrium length 
of the spring as indicated for three different values of b in 
Fig. [3] Consequently the critical value u C 2, which marks 
the transition to finite values of /3, is also sensitive to the 
bead separation. This monotone ^-dependence of u C 2 is 
shown for three different angles a in Fig. |4j 
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FIG. 4: The critical velocity u C 2 for the onset of the stair-like 
motion is shown as function of 6 for three different inclination 
angles a and u = 1. 

During its stair-like motion in the xy plane the dumb- 
bell changes its orientation, which is described by the 
angle 5 between r 12 and the x axis, while passing the 
saddles of the potential as indicated by the lower panel 
in Fig. [5] The amplitude of 6 depends on the dumbbell 
length b. With increasing values of u and a those reori- 
entation steps occur more often, accordingly the typical 
frequency of 8{t) increases as well as the angle [3. 

Alternatively one may vary the inclination angle a in- 
stead of the flow velocity u. If u is fixed at values larger 
than Ud, the drag force is too weak at a small an- 
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FIG. 5: The trajectories of the two beads in the xy plane 
are shown in the upper panel for the parameters a — 0.2, u — 
1, b = 1.5. The lower panel shows the corresponding evolution 
of the orientation angle 5 of the dumbbell. 

gle a to move the dumbbell across the potential barriers 
in the y direction. But above a critical value a c , which 
again depends on the length b of the dumbbell, one finds 
the transition from the motion along a potential valley 
to stair-like trajectories. The relation between (3 and the 
inclination angle a is shown in Fig. [6] Note, if the HI is 
neglected, all curves (3(a) for dumbbells of different sizes 
coincide with the one of a single bead (HI, which means 
that without HI the deflection angle has no ^-dependence. 
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FIG. 6: The deflection angle (3 is shown as a function of the 
flow inclination angle a for different values of the dumbbell- 
extension (dotted: 6 = 3, dashed: 6 = 1, dash-dotted: 6 = 
0.5) and u = 1. The solid line is obtained for a dumbbell in 
the free-draining limit (without HI) and does not change with 
6. It coincides with the curve /3(a) of a single bead (crosses). 
The inset shows the critical angle a c versus 6 for a dumbbell 
with HI (solid) and without HI (dashed). 

Just below the critical values a c and u C 2 thermal fluc- 
tuations support the passage of the beads between neigh- 
boring potential trenches and therefore may induce a 
stair-like motion in this parameter range. We investi- 
gated this influence of the thermal noise by taking into 
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FIG. 7: The deflection angle j3, averaged over 500 independent 
runs, of a dumbbell with b = 1 as a function of the flow angle 
a at u = 1 for different noise amplitudes. 

account the additive stochastic force Ff in the equation 
of motion ([I]). 

The deflection angle (3 as a function of a is shown for 
different values of the noise amplitude in Fig. [7J These 
numerical results, which are obtained by averaging over 
500 independent runs, show that the dumbbell motion is 
most sensitive to thermal fluctuations close to the critical 
angle a c or close to the critical velocity u C 2- The onset 
of the stair-like motion is shifted to smaller values of u 
or a as explained above, and in the limit of small noise 
amplitudes the deterministic result is approached. In the 
case of a noise amplitude in the order of Vq, the angles a 
and P become similar, indicating that the thermal motion 
is strong enough to kick the beads across the saddles 
of the potential. In this case all dumbbells are equally 
deflected and particle sorting is no longer possible. 

Another interesting dependence of the deflection angle 
on the noise strength is shown in Fig.[8l where (3 exhibits 
an interesting resonance like behavior as a function of the 
bead distance b. The deflection angle is clearly reduced 
in the regime, where the dumbbell length is a multiple 
of the wavelength and a higher excitation energy of the 
dumbbell is required. According to the behavior shown 
in Fig. [8l dumbbell sorting with respect to the size is 
most efficient within the range 0.5 < b/X < 1.5 in the 
limit of k B T < Vq. 

Conclusions and discussion.- We investigated the flow 
induced transport of mesoscopic and dcformable dumb- 
bells through a spatially periodic potential and found 
three different regimes of motion. The transitions be- 
tween these regimes depend sensitively on the flow veloc- 
ity, on the size of the dumbbell, and therefore on the hy- 
drodynamic interactions between the two beads. Below 



a first critical velocity the dumbbell is locked by the po- 
tential while above this threshold the two beads start to 
move along a potential valley. The critical flow strength 
depends on the length of the dumbbell in a periodic way. 
If the flow velocity increases further at a finite inclination 
angle and crosses a second threshold, then a transition 
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FIG. 8: The deflection angle /3, averaged over 500 independent 
runs, shown as function of the dumbbell size b. It exhibits an 
interesting resonance like behavior, especially for the noise 
amplitude fcsT/Vb = 0.1. Parameters: u = 1, a = 0.2. 

to a stair-like motion occurs. The preferred mean direc- 
tion of this motion is neither along the flow nor along 
a symmetry direction of the potential, and the resulting 
deflection angle depends on the length of the dumbbell. 
This allows for the use of periodic potentials, with their 
symmetry axis oblique to the flow direction, to deflect 
particles with respect to their size. The strongest sen- 
sitivity of the deflection angle on the dumbbell length, 
which may be employed for particle sorting purposes, oc- 
curs in the cross over regime to the stair-like motion. This 
cross over is sensitive to thermal noise because thermal 
fluctuations make the transition less sharp. Moreover 
they induce a resonance like structure in the deflection 
angle caused by the interplay between the wavelength of 
the potential and the size of the dumbbell. These ef- 
fects, which weaken the cfnccncy of particle sorting, may 
be suppressed by increasing the potential amplitude and 
the flow velocity simultaneously. Note, that the results in 
the regime of the transition to the stair-like motion can 
be obtained qualitatively too if an experimentally eas- 
ier accessible potential is chosen by neglecting the first 
contribution on the right hand side of Eq. ([3]) . 

We would like to thank L. Holzer and S. Schreiber for 
instructive discussions. This work has been supported 
by the German Science Foundation through the priority 
program on micro- and nanofluidics SPP 1164. 



[1] P. T. Korda, M. B. Taylor, and D. G. Grier, Phys. Rev. [2] M. De Luca, A. Mazzoni, and A.Vulpiani, Europhys. 
Lett. 89, (2002) 128301 . Lett. 60, (2002) 181 . 



5 



[3] M. P. MacDonald, G. C. Spalding, and K. Dholakia, Na- 
ture 426, (2003) 421. 

[4] J. A. Davis and et. al., Proc. Natl Acad. Sci. 103, (2006) 
14779 . 

[5] R. L. Smith, G. C. Spalding, K. Dholakia, and M. P. 

MacDonald, J. Opt. A: Pure Appl. Opt. 9, (2007) 134 . 
[6] B. Kaoui et al, Phys. Rev. E 77, (2008) 021903 . 
[7] S. Kapishnikov, V. Kantsler, and V. Steinberg, J. Stat. 

Mech. (2006) 01012 . 
[8] M. P. MacDonald and et. al., Science 296, (2002) 1101 . 
[9] K. Ladavac, K. Kasza, and D. G. Grier, Phys. Rev. E 

70, (2004) 010901 (R) . 
[10] J. Bammert, S. Schreiber, and W. Zimmermann, Phys. 



Rev. E 77, (2008) 042102 . 
[11] P. M. Johnson, C. M. van Kats, and A. van Blaaderen, 

Langmuir 21, (2005) 11510 . 
[12] T. C. B. McLeish and R. G. Larson, J. Rheol. 42, (1998) 

81 . 

[13] J. Rotne and S. Prager, J. Chem. Phys. 50, (1969) 4831 

[14] J. K. G. Dhont, An Introduction to dynamics of colloids 

(Elsevier, Amsterdam, 1996) . 
[15] M. Pelton, K. Ladavac, and D. G. Grier, Phys. Rev. E 

70, (2004) 031108 . 



